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ABSTRACT 



Context. Thermally pulsating asymptotic giant branch (AGB) stars are the main producers of slow neutron capture (s-) process 
elements, but there are still large uncertainties associated with the formation of the main neutron source, '^C, and with the physics of 
these stars in general. Observations of .s-process element enhancements in stars can be used as constraints on theoretical models. 
Aims. For the first time we apply stellar population synthesis to the problem of .s-process nucleosynthesis in AGB stars, in order to 
derive constraints on free parameters describing the physics behind the third dredge-up and the properties of the neutron source. 
Methods. We utilize a rapid evolution and nucleosynthesis code to synthesize different populations of ^-enhanced stars, and compare 
them to their observational counterparts to find out for which values of the free parameters in the code the synthetic populations fit best 
to the observed populations. These free parameters are the amount of third dredge-up, the minimum core mass for third dredge-up, 
the effectiveness of '''C as a source of neutrons and the size in mass of the '^^C pocket. 

Results. We find that galactic disk objects are reproduced by a spread of a factor of two in the effectiveness of the '''C neutron source. 
Lower metallicity objects can be reproduced only by lowering by at least a factor of 3 the average value of the effectiveness of the 
'^C neutron source needed for the galactic disk objects. Using observations of i-process elements in post-AGB stars as constraints we 
find that dredge-up has to start at a lower core mass than predicted by current theoretical models, that it has to be substantial (A > 0.2) 
in stars with mass M < 1.5 M© and that the mass of the pocket must be about 1/40 that of the intershell region. 

Key words, stars: AGB and post-AGB - stars: abundances - nuclear reactions, nucleosynthesis, abundances 
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About half of the elements heavier than iron are produced by 
slow neutron captures (s-process) in the deep layers of thermally 
pulsating (TP-) asymptotic giant branch (AGB) stars, which rep- 
resent the late evolutionary phase of objects with initial mass 
, M < 8 Mq. These stars have degenerate carbon-oxygen cores 
■ surrounded by two burning shells, and undergo a cyclic in- 
' stability of the helium-burning shell known as thermal pulses. 
5 -process elements are believed to be synthesized in the thin 
helium-rich layer between the helium-burning and hydrogen- 
burning shells known as the inter-shell region, or intershell. 
After each thermal pulse the hydrogen-rich convective enve- 
lope may penetrate into the intershell and bring the synthesized 
heavy elements, together with the products of helium burning 
(mainly carbon), to the stellar surface in a phenomenon called 
third dredge-up (TDU). Strong stellar winds gradually ablate the 
envelope and enrich the interstellar medium with the nucleosyn- 
thesis products. Many of the physical processes taking place in 
TP-AGB stars are still highly uncertain, such as how the main 
source of free neutrons required for the .s-process is produced 
(see below), under which conditions TDU occurs and how effi- 
cient it is, and the strength and driving mechanism of mass loss. 
Many of these uncertainties are related to our poor understand- 
ing of mixing processe s in stars. For a recent review on AGB 
stars see iHerwid (l2005h and for a discussion on u ncertainties in 
AGB stellar models see, e.g.. lLugaro et al.l ( l2003h . 



The synthesis of .s-process elements requires a source of 
free neutrons that can be captured by iron seeds to build up 
heavier nuclei along the isotopic stability valley. During the 
period of quiescent hydrogen shell burning in between ther- 
mal pulses, the inter-pulse period, the conditions in the inter- 
shell (temperature, density and helium abundance) favour the 
^^C(a, n)"'0 reaction. For this reaction to be an effective source 
of free neutrons a '^C-rich region or pocket in the intershell 
is needed. The production of '^C takes place by means of 
the ^^C(p,yy^N(J3^vy^C reaction. It is assumed that the pro- 
tons needed for this reaction penetrate the top of the inter- 
shell after TDU has occurred, as a result of a mixing process 
of as yet uncertain nature: semi-con vection (jHo llowell 
Il988 ), hydrodynamical overshooting ( Herwig et alj|1997i) . rota- 
tion dLanger et al.ll 19991; iHerwig et al] |2003!) and gravity waves 
dDenissenkov & Tout! 2003h have all been proposed. As a con- 
sequence large uncertainties exist in the size in mass of the '^C 
pocket and the '^C abundance profile inside it. In addition, due 
to the existence of several neutron poisons, especially '"*N, not 
all free neutrons released by the '^C(a, n)'^0 reaction participate 
in the i-process, so the effectiveness of ^^C as a neutron source 
for the i-process is also very uncertain. Another neutron burst 
is also released in the convective intershell during the thermal 
pulses due to the '^^Ne(a, n)^^Mg reaction, but the activation of 
this reaction is only marginal, which makes its contribution to 
the neutr on exposure small er than that of the ^^C(a, n)'^0 reac- 
tion (see Bus so et aLlll999l for a review of i-process nucleosyn- 
thesis in AGB stars). 
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Among stars that are enhanced in .s-process elements (s- 
enhanced stars) one distinguishes intrinsic and extrinsic objects. 
Intrinsic i-enhanced stars are typically late-type giants of spec- 
tral class S and C that show over-abundances of .^-process ele- 
ments (see, e.g.. ISmith & Lambert! 1 19901: Abia et al. in- 
cluding the radioactive element Tc (lMerrilllll952r With a half- 
life of a few million years, the presence of Tc lines indicates 
that the synthesis of i-process elements happened recently and in 
situ. Therefore intrinsic .s-enhanced stars must be objects in the 
TP-AGB phase or post-AGB phase, which have recently expe- 
rienced thermal pulses and TDU. The study of elemental abun- 
dances in these stars gives clues on the uncertain physical phe- 
nomena that take place in their interiors. Extrinsic .s-enhanced 
stars also show ove r-abundances of 5-pr ocess elements, but no 
Tc lines are present dJorissen et al.ll 19931) . This indicates that Tc 
has decayed since the s-process elements were synthesized. In 
addition, extrinsic stars are often observed to be in an evolution 
phase earlier than the TP-AGB, i.e., they are (sub)giants or main 
sequence stars. This implies that the s-process elements were not 
produced in situ, but that they have been accreted from a more 
massive TP-AGB companion star. These stars can therefore act 
as probes for studying mass transfer processes in binaries and 
for tracing the nucleosynthesis that occurred in their companion 
starS; 

iBusso et"aLl (12001 1) carried out a detailed comparison be- 
tween observations of s-process enhanced stars and model pre- 
dictions for the s-process in single AGB stars of different metal- 
licities. They found that a large spread (a factor of ~ 20) in the 
effectiveness (see ^23] for a definition of this parameter) of the 
main neutron source '^C was needed to match spectroscopic ob- 
servations at a given metallicity. Their analysis is based on a 
small set of initial stellar masses and they have taken into con- 
sideration only the final abundances of their stars. In this pa- 
per we improve on their analysis by comparing observations of 
different types of stars belonging to the AGB family to stellar 
population synthesis models computed with the inclusion of the 
s-process. Our use of a rapid synthetic evolution code allows us 
to study a large set of finely spaced initial masses and metallici- 
ties and to trace the complete AGB evolution of each star. Using 
this method we are able to put much tighter constraints on the ef- 
fectiveness of '^C as a neutron source. In addition, by comparing 
our models to post-AGB observations we also put constraints on 
the minimum core-mass for TDU, the TDU efficiency and the 
size of the '-^C pocket. 

The AGB stellar models and the associated free parameters 
are described in ^ Our stellar population synthesis method is 
described in ^ and the results are presented in ^ In ^ we 
discuss our results and draw our conclusions. 



2. The AGB models 

Our TP-AGB star models are calcula ted with a modified ver- 
sion of the r apid evolutionary code by iHurlev et al.l (1200 0) and 
llzzard et"aLl (12004) . We limit ourselves to a brief overview of the 
ingredients of this code and a discussion of the most important 
free parameters. A detailed description of our modifications and 
improvements can be found in Appendix|A] Stellar evolution and 
nucleosynthesis are modelled synthetically, i.e., b y means of an- 
alytica l fits to the detailed evolutio nary models of iKarakas et al.l 
(I2002h and IStancliffe et akl (|2004 . and in some cases by tabu- 
lar interpolation of these detailed models. Fits and interpolations 
depend on global stellar parameters, such as metallicity, mass 
and core-mass. This synthetic approach makes the calculation of 
the evolution very fast. The s-process element nucleosynthesis is 



calculated by interp olating results based on detailed models by 
iGallino et al.l ( Il998l) . as described below. With these models we 
are able to follow the chemical abundances of stars as a func- 
tion of time, from their initial abundances at the beginning of the 
main s equence, for which we use those of lAnders & Grevessd 
(1 19891) . to their chemically enhanced abundances of the latest 
evolutionary TP-AGB and post-AGB phases. 

Many uncertainties exist in the detailed evolutionary models, 
which we treat as free parameters. In the following subsections 
we discuss the free parameters that are most relevant for stellar 
chemical evolution. 



2. 1 . Minimum core-mass for tinird dredge-up 

Detailed models of TP-AGB stellar evolution which find 
TDU show that it only occurs when the stellar core-rnass is 
greater than about 0.6 Mq (e.g ., Boothrovd & Sackmannlll988l : 
lLattanziolll989HStraniero et al.lll997h . but there is disagreement 
of a bout 10% on the exa ct value of this minimum core-mass (see, 
e.g..'Lu garo et al.l2003l) . The luminosity and surface abundances 
of carbon and s-process elements in TP-AGB stars of mass M < 
2 Mq depend strongly on this minimum core- mass. We have 
model led th e TDU by using the fit to the models of lKarakas et al.l 
(l2002h from llzzard et al.l (|2004|) . which includes the free param- 
eter AM™". This parameter offsets the minimum c ore-mass by 
a fixed amount compared to the results of Ka rakas et al.l (|2002|) . 
allowing stars of low mass to undergo TDU even if they do not 
do so it in the detailed models. 



2.2. Tinird dredge-up efficiency 

Detailed evolution models disagree on the amount of third 
dredge-up that occurs in a star of a particular mass and metal- 
licity, i.e., how d eep the convective envel ope penetrates the in- 
tershell (see, e.g.. lFrost & Lattanzioll 9961 and rMowlavill 19991 for 
discussion). The amount of dredge-up can be measured by the 
third dredge-up efficiency 



A = 



AMtdu 
AMh ' 



(1) 



where AMjou is the reduction in mass of the hydrogen ex- 
hausted core as a result of dredge-up and AMh is the amount 
of mass by which the hydrogen-exhausted core increased due to 
hydrogen burning during the previous inter-pulse period. Thus, 
if A - there is no dredge-up at all, while if A - I there 
is no net growth of the hydrogen-e x hauste d core. We apply 
the parametrization from llzzard et"al] (2004') in which A grows 
asymptotically with pulse number to a maximum value 



(2) 



Here A^^^ is a function of mass a nd metallicity which fits the 
^max value of the detailed models d Karakas et al.ll20()2|) and the 
free parameter /l„iin is used to set a lower limit. In this way stars 
which show negligible TDU efficiency in the detailed models are 
allowed to have a substantial TDU efficiency when Amin is set to 
be larger than zero. 

2.3. Effectiveness of ^^C as a neutron source 

At the end of a TDU episode protons from the bottom of the 
H-rich envelope may penetrate beyond the chemical boundary 
into the helium- and carbon-rich intershell. They react with '^C 
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to produce '^N, which decays into '^C forming the '^C pocket. 
If there are enough protons, '^C will also capture protons and 
produce ''*N. The '^C in the pocket is a source of neutrons by 
means of the ^^C{a,ny^O reaction, but it has to compete with 
'"^N as a sink of free neutrons due to the ''*N(n, reaction, 
which has a relatively high neutron capture cross section. As a 
result of this competition an effective local '■'C abundance can 
be defined, which parametrizes the number of free neutrons per 
heavy seed nucleus in the intershell available for the .s-process 
to take place. The abundance ratios of i-process elements in the 
'^C pocket are determined by this local '^C abundance, which is 
scaled relative to a standard case. This standard case was defined 
by Gallino et al. (1998) and Busso_etal,) (|2001|), based on the '^C 
pocket found by|Hpllowell & Iben ( 1988) in their detailed evolu- 
tionary calculations. G allino et alj (il998<) showed that with this 
pocket low mass AGB stars of half solar metallicity produce the 
main i-process component for the Sun. We denote this scaled 
effective '''C abundance as '^'Ceff, so that the standard case cor- 
responds to '■'Ceft' = 1. 

Comparing the absolute elemental abundances produced by 
s-process models directly to observations is difficult because 
they depend on the dilution of material from the intershell into 
the envelope, which in turn depends on several uncertain factors 
such as the size in mass of the '^^C pocket, the amount of TDU, 
mass loss and mass accretion. However, the i-process element 
abundance ratios remain mostly unaffected by these processes, 
so they provide direct constraints on '"'Ceff (see 



2.4. Size of the "C pocket 

In order to calculate the abundances in the stellar envelope we 
must know the element mass fraction in the intershell at the mo- 
ment when the star undergoes TDU. For i-process elements, we 
calculate these inte rshell abundances by using the results of de- 
tailed models from lGalUno et al.l (Il998l) . which consist of a grid 
of intershell elemental abundances at the time of TDU for stel- 
lar masses of 1.5, 3 and 5 Mq, metallicities lO""* < Z < 0.02, 
effective local '■'C abundance (1/24 < '■'Ceff < 2) and number 
of thermal pulse followed by TDU (2 < A^tdu < 29). We know 
the mass of the '^C pocket and of the intershell in these detailed 
models, which we use to trace back the local element mass frac- 
tions in the '^C pocket just before the thermal pulse occurs. We 
map a grid of .s-process element mass fractions in the '^C pocket 
for different stellar masses, metallicities, '^Ceif and pulse num- 
ber, which is independent of the choice of the size of the '''C 
pocket and that of the intershell. We interpolate linearly on this 
grid and if the stellar mass and/or pulse number needed for the 
calculation is out of the grid bounds we use the result of the 
closest point on the grid. 

Note that with this procedure we do not properly account for 
the ^^Ne neutron source. However, this neutron source does not 
contribute to defining the overall i-process distribution, at least 
for stellar models of mass M < 4 Mq, which are believed to be 
the counterparts of the observed .^-process enhanced stars that 
we will discuss in ^(e.g., Abia et a l. 200]^. 

We assume that the '^C pocket mass is a fixed fraction, 
/i3c of the mass of the intershell at the moment just before 
TDU occurs, Mjs, which in turn is well approximated by the 
maximum mass of the convective intershell pocke t dur i ng the 
thermal pulse. We employ the relation proposed bv llbeni (1 19771) 
for the maximum convective intershell mass. 



where is the mass of the H-exhausted core. This relation 
is derived for core masses larger than O.95M0, but when ex- 
trapolated to lower core masse s it yields sim i lar in tershell mass 
values to those considered by iGallino et al.l (Il998l) in their de- 
tailed models. The value of /"cis can be constrained by obser- 
vations of surface abundances of i-process elements, although 
these abundances also depend on the amount of dredge-up that 
the star has experienced. This degeneracy is broken by studying 
the abundances of other elements that are not synthesized via the 
i-process, as discussed in 94.1.41 

2.5. Mass loss 

Mass loss is an uncertain factor in AGB evolution which affects 
the enhancement of i-process elements on the surface of a TP- 
AGB star (Straniero et al. 2003). The mass loss rate determines 
the stellar lifetime as a TP- AGB star by truncating i t when the 
envelope is lost. In the models of iGallino et alj (Il998i) the inter- 
shell 5-process abundance ratios vary from one pulse to another, 
until after roughly 20 pulses they reach an asymptotic value. The 
asymptotic values are usually considered to be the i-element ra- 
tios of a t ypical pulse and are often compared to the observed ra- 
tios (e .g., OBusso et al .112001 1: 1 Straniero et al.ll2003l:lRevniers et al.1 
|2004|) . Depending on the mass loss rate applied in the models, a 
star with mass M < 1 .5Mq can undergo from very few thermal 
pulses (or none) to several tens, which affects its surface compo- 
sition and yield during the T P- AGB phase. In this work we use 
the mass-loss prescription bv lVassiliadis & WoodI (Il993h : 



log 



10 



Mo yr-' 



-11.4+0.025 



P 

day 



100 max 



M - 2.5Mo 



(4) 

where log P/day = -2.07 -i- 1.94 log R/Rq - 0.9 log M/Mq is 
the logarithm of the Mira pulsation period. We do not attempt to 
vary the AGB mass-loss rate, but we discuss the consequences 
of this assumption on our results in ^ 

The mass-loss rate during the giant branch phase affects the 
envelope mass at the moment when a star begins its TP-AGB 
phase. We have modelled mass- loss from giant branch stars with 
the widely used prescription bv lReimersI ( 1975h . 



Mgb = 4 X 10-'3/,GB (i/Lo) (R/Ro) (M/Mo)-' yr"', (5) 

with //GB =0.3, which is within the limits set by the morphol- 
ogy of horizontal branch stars in the gala ctic globular clusters 
(llben & Renzinilll983l;ICarraro et al.lll996l) . 



3. Population synthesis 

We make use of our synthetic evolutionary models to carry out 
population synthesis of single stars, on a grid of 50 logarithmi- 
cally spaced metallicities ranging between 10 < Z < 0.025 
and 500 logarithmically spaced initial masses, whose range de- 
pends on the metallicity and on whether we are modelling a pop- 
ulation of intrinsic or extrinsic i-enhanced stars. 

In the case of intrinsic stars, for each metallicity we apply 
the initial mass ranges shown in Fig.[T] with 8 M© as a maxi- 
mum. We derive these ranges from combining our stellar evo- 
lution models with the age -metal licity relation (AMR) for the 
galactic disk of iPont & Eve^ (|2004 



logio(Mis) = -1.835 + 1.73Mc - 2.67MI 



(3) 



Fe 

h" 



logi 



f- 17.82 
-23.24 



: 0.25 dex. 



(6) 
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-0.8 -0.6 -0.4 -0.2 
[Fe/H] 



0.2 



Fig. 1. Range of initial masses, Mi, of our population synthesis 
grid as a function of metallicity, depicted by the shaded area. 
For a given metallicity this shows the initial-mass range of stars 
which can be observed in the TP- AGB phase, with a maximum 
mass of 8 Mq. 



where the notation for the element number density ratio 
[y] = logio(^/^j is used and f is the age in Gyr, with the 
age of the Universe (13.7 Gyr) as its upper limit. Then, from our 
single-star models we consider the range of initial masses for 
which a star spends at least some of its TP-AGB life within the 
age range allowed by the AMR. The probability for an intrinsic 
star to have a given abundance is weighted by the time it spends 
showing this abundance and by its contribution ac cording to the 
initial mass function (IMF). We use the IMF, ^, of iKroupa et al.l 
(Il993h . given by 



0.035Mri-3 
0.019Mr2-2 
0.019Mr2-7 



if 0.08 < Mi <0.5, 
if 0.5 < Mi < 1.0, 
if 1.0 < Mi, 



(7) 



where Mi is the initial stellar mass in solar units. 

Given that we employ single stellar models, we treat the case 
of extrinsic stars indirectly. These stars have acquired their s- 
process enhancements by accreting part of the mass ejected by 
an initially more massive binary companion. Therefore we as- 
sume a simple model in which the initial mass of the companion 
is assumed to be distributed according to the single-star IMF. 
Hence the probability of an extrinsic star to show a given ele- 
mental abundance is weighted by the amount of mass which is 
lost by its companion in the form of that element and by the IMF 
weight of this companion. Because mass transfer could have oc- 
curred at any time since the formation of the system, we consider 
all stars whose evolution is beyond the TP-AGB phase, which in 
each metallicity bin corresponds to initial masses greater than 
the lower limit in Fig. [T] and smaller than 8 Mq. For the mod- 
els with [Fe/H] < -1 we consider 0.8 M© to be the lower initial 
mass limit. 



4. Results 

In this section we present the results of our population synthesis 
models and compare them to different sets of observational data. 
The observations provide measurements of the surface abun- 
dances of .s-process elements relative to that of iron. Following 
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Fig. 2. Envelope [hs/ls] during the evolution of 1.5 Mq and 3 
Mq star models of different metallicity, as indicated in the panel, 
and with '^Ceif=l. 



the convention, we calculate for each star in our grid its heavy 
i-process (hs) element abundance ratio. 
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and light i-process (Is) element abundance ratio. 



ii 
Fe 



Y 

Fe 



Zr 

Fe 



(8) 



(9) 



These ratios are altered by dilution in the stellar envelope, which 
is affected by several uncertain factors such as the size in mass of 
the '^C pocket, the amount of TDU, mass loss and, in the case of 
binaries, mass accretion. We also study the ratio of hs-elements 
to Is-elements, 



hs 




hs 




Is 


Is" 




Fe 




Fe 



(10) 



This ratio reaches asymptotically the [hs/ls] values of the inter- 
shell after a number of pulses, making it mostly unaffected by 
dilution. 

The asymptotic envelope ratio is therefore only determined 
by ^^Ceff (see ^2.3b . but not by other free parameters in our 
model, such as the size of the '^C pocket {fac.is)- Fig- El 
shows enve l ope [h s/ls] ratios as obtained from the models of 
iBusso et" ar ( 2001 ), which in each model reach their correspond- 
ing intershell asymptotic value, but take a number of pulses to 
do so. Hence, considering stars that have experienced only few 
TDU episodes, either because they are in an early stage of the 
TP-AGB phase or because they have masses M < 1 .5 Mq, re- 
sults in a natural spread of [hs/ls] values. Only in AGB stars 
with M > 5 Mq is the amount of dredged-up material is so small 
that the surface [hs/ls] ratios differ from those of the intershell. 

To compare these abundance ratios with the observed data 
we have binned the resulting number distribution with a resolu- 
tion of 0. 1 dex. Whenever a star has a certain abundance ratio, its 
probability is added to the corresponding bin. This is done sepa- 
rately for the different metallicity bins and in each we normalize 
the resulting probability distribution by dividing it by its highest 
value. 

We select to be i-enhanced those stars which have [Is/Fe] 
and/or [hs/Fe] larger than 0.1 dex. First we analyze the [hs/ls] 
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Intrinsic MS/S stars (CI 3g„o=4/3) Intrinsic MS/S stars (CI 3effo=2) 




-1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 

[Fe/H] [Fe/H] 



Fig. 3. Intrinsic MS/S stellar population synthesis results compared to the observations. The panels show our results calculated 
with four increasing '^Cefr values as indicated above the panels. The grey scale is a logarithmic meas ure of the normalized number 
distribution of stars over [hs/ls]. The crosses are the observational data gathered bv lBusso et alj ( 1200 ll) . which have an average error 
given by the size of the cross in the upper right of each plot. 



ratios of intrinsic i-enhanced stars, which have a strong depen- 
dence only on the '■'Cefr value. In all the runs we use AM™'" = 
-O.O65M0, /Imin - 0.2 and fnc.is - 1/40, based on our attempts 
to fit the post- AGB star observations of [Zr/Fe] ratios and carbon 
abundances, as explained in detail in 94.1.41 In 94.2l we study the 
[hs/ls] ratios, but this time of extrinsic .s-enhanced stars, which 
give clues on ' ''CefF at metallicities lower than that of the galactic 
disk. 



4.1. Intrinsic s-enhanced stars 
4.1 .1 . [hs/ls] ratios of MS/S stars 

MS/S stars are late-type giants that show lines of ZrO and TiO, 
indicating that they are enhanced in s-element abundance, have 
low surface temperature and their C/O ratio is smaller than unity. 
To account for these stars we select from our population synthe- 
sis i-enhanced stars in the TP-AGB phase of evolution with a 
surface effective temperature less than 3500 K and C/O < 1 in 

ratio to ob- 



Busso et aTl 



their envelope. We compare our results for the [hs/ls 

servati onal data of intrinsic MS/S stars gathered by[ 

(1200 lb . In Fig. |3] our results only appear at [Fe/H] > -0.6 be- 
cause low-metallicity stars that under-go TDU episodes are too 
hot to be classified as MS/S objects. This is consistent with 



the observations. Increasing '^Ceff in our models increases the 
[hs/ls] ratios on the low-metallicity side of the plot, but at high 
metallicities the [hs/ls] ratios are decreased. This is explained 
as follows: If the number of heavy iron seeds is decreased for 
a given number of available free neutrons, then the [Is/Fe] ra- 
tio increases due to i-process reactions and [hs/ls] decreases, as 
is observed at high metallicities in our results. But if [Fe/H] is 
decreased even more then the Is-element synthesis saturates and 
the hs-element abundance increases with respect to iron, rais- 
ing both ratios [hs/Fe] and [hs/ls], as can be seen in our mod- 
els at low metallicities (Travaglio et al. 1999; Busso et al. 200 1'). 
Thus '■'CefF basically determines at which metallicity our mod- 
els will show [hs/ls] ratios below solar and where they will rise 
to solar and beyond. This effect cannot be seen clearly in our 
MS/S star synthesis due to the lack of .s-enhanced MS/S stars be- 
low [Fe/H] 0.6, but is common to all the stellar populations 

that we have synthesized and is more evident in SC and C stars 
and post- AGB stars ( g4.1.2| and 94.1.31 respectively). In contrast 
with other studies where only the asymptotic values of [hs/ls] a re 
considered (e.g.. lGorielv & M owlavi"2000!: lBusso et al.ll2001h . a 
natural spread in this ratio is observed in our results. This is ex- 
plained by our tracing of the whole TP-AGB evolution, where 
the initial dozen pulses have different i-process ratios. Stars of 
mass M < 1.5 Mq, which are the most numerous due to the IMF 
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error given by the size of the cross in the upper right of each plot. 



and which do not undergo more than 10-20 pulses, contribute 
most to this effect. Fig. [3] shows that the best fit to the observed 
distribution of MS/S stars is obtained with ^"'Ceff ~ 4/3 and a 
spread in '^Ceff values is not needed to match all the observa- 
tions. Given that the error bars on the individual observations are 
large (as indicated in Fig. [3]), and that the observed [hs/ls] values 
are close to 0, the value of '^Ceff is not strongly constrained by 
these observations. The (on average) negative values of [hs/ls] 
are matched, within observational errors, by the models if '^Ceff 
is between 2/3 and 2. As the [hs/ls] ratio is not sensitive to di- 
lution, changes in the other free parameters do not affect these 
results. 



4.1 .2. [hs/ls] ratios of SC and stars 

SC and C stars are late-type giants that show lines of carbide 
molecules in their spectrum, which indicates that there is more 
C than O by number In the case of the SC stars, the C/O ratio is 
close to unity, given that they also show ZrO lines. To reproduce 
these stars we select from our population synthesis results those 
s-enhanced TP- AGB stars which have C/O > 1. We compare 
our synt hetic [hs/ls] r atios with t he observa t ional data gathered 
by Buss o et al.l (1200 ll) and from lAbia et alj (|2002|) as shown in 
Fig. |4] The spread in the [hs/ls] ratio is explained in the same 



way as for MS/S stars. The spread is smaller at low metallici- 
ties due to the small range of initial stellar masses considered 
at these metallicities resulting from the AMR (see Fig. [1]). The 
effect of changing '^Cefr on our synthetic [hs/ls] ratios is clearer 
in this population of stars than in the MS/S stars, particularly for 
[Fe/H] < -0.4. As in the case of the MS/S stars, the other free 
parameters do not affect the [hs/ls] ratios. Most of the observed 
SC and C stars cluster around [Fe/H] ^ and slightly nega- 
tive values of [hs/ls]. These stars can be fitted with any value of 
'^Ceff. The diagonal crosses in Fig.|4]are halo C stars suspected 
to be intrinsic. They are not expected to follow the same AMR 
as those in the galactic disk, but the fact that they are old and 
of low mass is approximately modelled by our AMR. Three of 
these stars are matched by our models for values of ^^C^e > 2/3, 
but a lower '"'Ceff cannot reproduce their [hs/ls]. The lack of stars 

below [Fe/H] 1 in our models is due to the fact that, with our 

adopted AMR, lower metallicity stars do not become TP-AGB 
stars within the age of the Universe. This effect is observed in 
all our results involving intrinsic stars and is consistent with the 
observations. 
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4.1 .3. [hs/ls] ratios of post-AGB stars 



Post-AGB stars are in the fast evolutionary phase between the 
AGB and white dwarf track. They suffer from strong mass loss 
and are hot, but not hot enough to ionize their circumstellar 
medium. We select them from our models by choosing those s- 
enhanced TP-AGB stars that have less than O.OSMq left in their 
envelopes. We compar e our synthet i c pop u lation to the observa - 
tional data gathered bylBusso et al.1 (llOOlh.lvan Winckell (l2003l) . 



tionai aata gatnerea bymusso et ai.l(izuuii) . ivan wincKeil (izuu 
iRevniers etall (|2004|) and iGiridhar & Arellano Ferrol (l2005h . 
Fig.|5]shows a spread, mostly at high metallicities, which is due 
to the different initial stellar masses taken into account in our 
population synthesis, the range of which becomes narrower to- 
wards low metallicities as we follow the AMR. Stars with M < 
1 .5 Mq experience a small number of TDU episodes. After their 
last thermal pulse (when they are observed as post-AGB stars) 
they have not yet reached the asymptotic [hs/ls] ratio, which con- 
tributes to the spread observed. Most of the data can be fitted 
with 2/3 < '^Ceif < 4/3, within the observational errors. There 
is only one exception at [Fe/H] ~ - 1 , which apparently needs a 
somewhat smaller '^Ceff. This object, IRAS07134-H1005, has a 
very high heliocentric velocity (Van Winckel & Reyniers 2000) 
suggesting that it belongs to the galactic halo rather than the 
galactic disk. 



4.1 .4. Zr and C abundances of post-AGB stars 

Post-AGB stars can also be used to provide constraints on the 
minimum core-mass for TDU, the TDU efficiency and the size of 
the '^C pocket. The observed values of [Zr/Fe] show a bimodal 
distribution: [Zr/Fe] ^ and 1 < [Zr/Fe] < 2, which suggests 
that some post-AGB stars did not experience TDU, while oth- 
ers di d and suffered a strong enhancement of Zr (van WinckeJ 
l2003h . This can be explained when we consider low mass stars 
(M < 1.5 M0) at sub-solar metallicities. They have low-mass 
envelopes which easily become strongly .s-enhanced. As shown 
in Fig. |6] only one TDU episode after the intershell material 
has become i-enriched is enough to raise their envelope abun- 
dance ratios by almost 1 dex. We compare the observa tional 
[Zr/Fe] ratios of Ivan Winckell (l2003h . IRevniers et all (12004 ) and 
IGiridhar & Arellano Ferrol (l2005b to our models, from which 
we consider all post-AGB stars, both i-enhanced and not s- 
enhanced. 

Intrinsically s-enhanced post-AGB st ars are observed a t 
metallicities as low as [Fe/H] ~ -1.0 dvan Winckell l2003h . 
When taking into account the AMR, this implies that at metal- 
licities [Fe/H] < -0.7 stars with initial masses above 0.9 Mq 
must experience TDU. At [Fe/H] < -0.3, the observed num- 
ber of .s-enhanced stars (which experienced TDU) is compara- 
ble to that of stars with [Zr/Fe] » (which did not undergo 
TDU). Thus, taking into account the IMF and the initial mass 
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range given by the AMR at low metallicities, we estimate the 
minimum initial mass at which TDU must take place so that the 
number of .s-enhanced stars is similar to that of non-enhanced 
stars. Applying this in our models, we find the correspond- 
ing minimum core mass for TDU as a function of metallicity. 
The resulting relation between minimum core mas s and metal- 
licity is best modelled by the fitting formula of .Karakas et alj 
(I2OO2I) . with an offset AM™" = -O.O65M0. This outcome de- 
pends somewhat on the adopted age-metallicity relation which 
sets the initial mass range of TP-AGB stars at each metallicity 
(Fig. [D> but the sensitivity of the mass of a low-mass star to its 
lifetime is only slight. Our result depends more strongly on the 



choice of mass loss. To test this we apply di fferent mass-los s 
rates on the TP-AGB with the prescription of ^ Reimersl ( Il975h . 
Eq. dU, by varying the coefficient jjcb from 0.5 to 3.0, and find 
that 0.055 < A M™" < 0.075. This result i s cons i stent with the 
finding s of, e.g., Groenewegen & de Jong ( 19931) . iMarigo et all 
(|1999[) . Ilzzard et all (|2004 and IStancliffe et all (|2005|) when 
studying the carbon star luminosity functions of the Small and 
Large Magellanic Cloud. 

We now address the question how much intershell material 
is dredged up to the surface. The surface i-process enhancement 
in our models is controlled by the dredge-up efficiency {A) and 
by the fractional mass of the '^C pocket relative to the inter- 
shell (/i3c,is)- The carbon enrichment in the envelope depends 
only on the amount of dredge-up, independent of the size of the 
'^C pocket. The observations show that below [Fe/H] ~ -0.3 
all i-enhanced post-AGB stars have C/O > 1 and that the non- 
enhanced ones have C/O < 1 dvan Winckell 120031) . We use this 
information to break the degeneracy between A and /ncis- Let 
Ns,noC be the observed number of .s-enriched post-AGB stars 
that have C/O < 1 and A^j be total number of .s-enriched post- 
AGB objects. Fig. I2] illustrates the ratio Ns,aoc/Ns as a function 
of metallicity calculated with our models for different values of 
Amin- With a simple application of Bayesian statistics we calcu- 
late that, given the observed data distribution, Ns,noc/Ns - and 
Ns - 11, there is a 32% probability that Ns,noc/Ns > 0.09 and 
only a 5% probability that Ns,noc/Ns > 0.23. The latter sets the 
95% confidence limit A^in ^ 0.2, as shown in Fig. |7] With this 
lower limit we calibrate /ucjs to fit the Zr enhancement of post- 
AGB stars. Fig.[8]shows our population synthesis results for the 
[Zr/Fe] ratio calculated with different values of A^i„ and /uc.is. 
compared to the observations. The upper left panel shows our 
results for A^in - 0. With this choice there are no i-enhanced 
post-AGB stars at low metallicities ([Fe/H] < -0.4) since these 
would all be low-mass stars due to the AMR and they experience 
negligible dredge-up. This is in contradiction to the observa- 
tions. Moreover, as shown above (Fig.|2ll, a choice of A„^in < 0.2 
is not compatible with the number of carbon enhanced stars. If 
we choose A^in = 0.2 (upper right panel) then low-metallicity s- 
enhanced post-AGB stars are indeed produced. The large spread 
in the i-enhancement at higher metallicities is due to stars with 
mass M > 3 Mq. However, the choi ce fup K = . 05, w hich 
roughly corresponds to that used by iGallino et al.l (Il998l) for 
their detailed nucleosynthesis calculations, gives too much Zr 
for many of the .s-enhanced post-AGB stars with [Fe/H] < -0.5. 
The lower right panel shows that most of the observations can 
be reproduced with A^i„ = 0.2 and /i3c,is ~ 1 /40, with two ex- 
ceptions which are extremely i-enhanced. From the constraint 
that all i-enhanced post-AGB stars are carbon stars, there is a 
68% confidence lower limit on A^i„ of 0.4. This implies that 
/"'CIS ~ 1 /lOO is needed to fit the post-AGB observations (lower 
right panel of FigO. However, this choice of /ncjs gives synthe- 
sized [Is/Fe] and [hs/Fe] ratios which are too small to reproduce 
those observed in MS/S, SC and C stars. Fig.|9]shows a compar- 
ison of the observed [Is/Fe] and [hs/Fe] of MS/S, SC and C stars 
to our models, calculated with AM™" = -O.O65M0, A^in - 0.2, 
/"CIS = 1/40 and '^Ces = 4/3. Most of the MS/S star abun- 
dances are fitted within the errors, with two high-metallicity ex- 
ceptions that are only matched with a larger choice of '^CeiF or 
/'c,is- The i-process abundances of most of the SC and C stars 
in the galactic disk are also matched well using the same choice 
of free parameters, except for three low metallicity objects that 
need a smaller '^Ceft . 

From the combined evidence of MS/S, SC, C and post-AGB 
stars, we find that with a few exceptions, all observations of in- 
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trinsic .s-enhanced stars can be matched by models with the fol- 
lowing set of parameters: 

- AM™° = -O.O65M0 

- Amin - 0.2 

- /'3c,is = 1/40 

- 2/3 < '^Ceif < 4/3. 

However, the constraints we have set here depend somewhat 
on our choice of other rather uncertain parameters, namely the 
mass-loss rate and age-metallicity relation (see §|5]l. 

4.2. Extrinsic s-enhanced stars 

Extrinsic i-enhanced stars can be in any evolutionary stage prior 
to the TP-AGB phase, but they are mainly giants and main se- 
quence stars. Their absolute i-process abundances cannot be 
studied without complete binary evolution models, but the abun- 
dance ratios of s-process elements are not substantially affected 
by binary processes. The [hs/ls] ratio represents that of the for- 
mer TP-AGB companion which produced the ^-process ele- 
ments and it can be reasonably approximated by taking into 
account the yields of single stellar evolution models. Fig. [TO] 
shows a comp arison of our synthe t ic [hs/lsl ra t ios to t he extrinsic 
star da ta from'Busso et al.' ('200l!). lAoki et all (l2002l) . lAbia et all 
(I2OO2I) and[^n Ecket al. (2003). Similar to our results for in- 
trinsic stars, the spread observed in these results is due to the 
range of initial masses that we use, although, in this case the 
mass range is wider and includes AGB stars up to 7-8 Mq. An 
AGB star of M > 5 M© has a very thin intershell and a massive 
envelope, so its surface [hs/ls] ratio varies smoothly between its 
value in the intershell and the solar ratio. 

We focus first on the galactic disk ([Fe/H] > -1) stars. Fig. 
[TOl shows that no single choice of '^Ceff value can make our re- 
sults fit all the observational data, but almost all observations can 
be fitted with a range of values 2/3 < '"'Ceff < 4/3. This is con- 
sistent with what we found for the intrinsic s-enhanced stars (see 

For [Fe/H] < -1, our synthetic [hs/ls] ratios are relatively 
insensitive to changes in '^Ceif because the hs-element synthesis 
saturates, opening the way to the synthesis of lead. There is some 
indication that a smaller value of '^Ceff is needed to reproduce 
objects with [Fe/H] < -2 which have a high [hs/ls], but even 
when using small '^Ceff values [hs/ls] does not exceed about 1 
dex. Thus we are unable to explain the extreme low-metallicity 
objects with [hs/ls] > 1.2. Our synthetic calculations do not ex- 
tend below [Fe/H] < -2.3 because of the lack of detailed models 
with metallicity below this value. 

4.2.1. Lead stars 

Lead stars have an enhanced lead abundance and are found at 
low metallicities. As discussed above, the [hs/ls] ratio remains 
roughly constant at low metallicities due to the saturation of the 
Is- and hs-el ement synthesis, wh ich indicates that lead is being 
synthesized (Gallino et alj|1998h . Consequently, the ratio of Pb 
to hs-elements is sensi tive to the choice of '^Ceff . We use [Pb/hs] 
data gathered by Biste rzo et al.l (|2006) to compare with our syn- 
thesis results. In the [Pb/hs] vs. [Fe/H] plane (Fig.[T2ll we see that 
our results show a pattern that shifts in metallicity if '^Ceff is var- 
ied, similar to the shift in [hs/ls] pattern at higher metallicities. 
The reason for this behaviour is the change in ratio of the num- 
ber of free neutrons to the number of seed nuclei, as explained 
in 94.1.11 We find that also for [Pb/hs] a spread naturally arises 
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Fig. 11. [Pb/hs] during the evolution of 1.5 Mq and 3 M0 mod- 
els of [Fe/H]=-2.3 with '^Ceff = 1/3. 



in our population synthesis r esults. This spread i s mostly due to 
the fact that, as discussed bv lGallino et al.l (11998 ) (see their Fig. 
6), the neutron exposure in the '^C pocket decreases with pulse 
number. This effect is even more pronounced at lower metallici- 
ties, causing the [Pb/hs] ratio to shift to much lower values as the 
ev olution proceeds. T his is illustrated by the models computed 
bv lBusso et al. l (l200Tl) . shown in Fig.[TT] As a consequence, there 
are two ridges of high probability in the model results shown in 
Fig. [12] One follows the higher [Pb/hs] values and corresponds 
to the contribution of low-mass stars (about 1 Mq), which is 
high due to their large IMF weight. The other ridge, which fol- 
lows intermediate values, is the contribution of stars with masses 
around 2.5 Mq. Despite their smaller IMF weight, the ejecta of 
these stars are massive enough to make their contribution com- 
parable to that of the low-mass stars. This effect is also visi- 
ble in our synthetic [hs/ls] results, e.g., in the top panels of Fig. 
[TOl Most observations are fitted by the contribution of ~ 2.5 
Mq Stars when choosing '^Ceif ~ 1/3 (lower left panel in Fig. 
[T2I 1. except for two outliers which have [Pb/hs] < 0. These out- 
liers, however, have l arge error bars which extend to [Pb/hs] x 
(Ivan Eck et al.ll2003h . The upper left panel in Fig. [T2]shows that 
a value of '^Cgfr as small as 1/12 still marginally fits the low- 
metallicity stars ([Fe/H] < -2). On the other hand, a choice of 
'■'Ceff ~ 2/3 results in too large [Pb/hs] values (lower right panel 
in Fig. [T2T l. This result reaffirms our results for the extrinsic- 
star [hs/ls] values which suggest that a somewhat smaller value 
of '^Ceff is needed to fit low-metallicity halo stars compared to 
those in the galactic disk. This is also consistent with what we 
found in ^4T3l regai'ding the object IRAS07134H-1005. 

As stated above our calculations do not extend below 
[Fe/H] < -2.3, but we notice that lead stars with lower metal- 
licity could also be explained if the trend shown by our mod- 
els continues at lower metallicities. This seems likely because 
[Pb/hs] does depend mostly on '^Cefi and not on other uncer- 
tain parameters in our models, however, explaining those stars 
at [Fe/H] < -2 which show [hs/ls] ^1.2 is still a problem. 

Finally, we note that ICui & Zhangj (j2006) also presented a 
possible solution to explain the spread of [Pb/hs] observed in 
low-metallicity stars. These authors obtained a spread of neutron 
exposures (and hence [Pb/hs]) in low-metallicity AGB stars due 
to variations of the inter-pulse period in stars of different masses. 
Their results are based on the assumption that the neutron irra- 
diation time is proportional to the inter-pulse period. However, 
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Fig. 8. Post-AGB stellar population synthesis results compared to the observations. The grey scale is a logarithmic measure of the 
normalized number distribution of stars in the [Zr/Fe] - [Fe/H] plane for '^Ceif = 2/3. As indicated above the panels, plots are 
presented for different values of and facis- The crosses are the observational data (see references in the text), which have an 
average error given by the size of the cross in the upper right of each plot. 



this assumption is incorrect since the neutron irradiation time 
depends instead on the timescale at which the '^ C(a, n)^^0 re- 
action s occur. At the temperature of 10** K that ICui & ZhangI 
(120061) use in their models, t he timescale for a ca ptures is of the 
order of 300 yr (see Fig. 2 of lGallino et al.|[T998 b. and this is in- 
dependent of the inter-pulse period. In our models, instead, the 
spread in [Pb/hs] derives naturally from the fact that the neutron 
exposure changes with time, especially at low metallicities. 

5. Conclusions and discussion 

Based on the results of our population synthesis study, we find 
that with the standard set of values for the free parameters in 
our AGB models it is not possible to reproduce the observed 
heavy-element abundances of i-enhanced stars. The minimum 
core mass for TDU has to be reduced by at least 10 per cent with 
respect to what theoretical models predict; low mass stars must 
have a significant TDU efficiency and the '^C pocket mass must 
be somewhat smaller than what is usually considered standard. 
We also find that the observations can be matched with a narrow 
range of '^'Ceff values and that the mean '^Ceif value in this range 
apparently decreases with metallicity. 

The observations in the metallicity range [Fe/H] > -1 are 
well fit by 2/3 < '^'Ceff < 4/3, i.e., a spread of a factor of 



2, while iBusso et all (l200lb needed a factor of ~ 20. This is a 
consequence of the natural spread that we find in the .s-process 
element ratios caused mainly by the truncation of AGB evolution 
due to the total loss of the stellar envelope. Stars with initial 
mass M < 1.5 Mq start the TP- AGB phase with a relatively 
similar core-mass, but different envelope mass. Consequently, a 
star with a certain low mass will experience fewer thermal pulses 
than one with a slightly larger mass, which affects their yields 
and composition during their life as a TP- AGB star. A different 
choice of mass loss in our models, either during the red giant 
branch or the AGB phase, would move the mass range in which 
this effect takes place, but not eliminate it. 

At [Fe/H] < - 1 a spread in [Pb/hs] arises because of the shift 
of the neutron exposure. Observations at these lower metallici- 
ties need a somewhat smaller value '^Ceff ~ 1/3, perhaps down 
to 1/12. This result may be due to the increasingly important 
effect of g-enhanc ed light neutron poisons at low metallicities 
jBusso et alj|1999h . which needs to be investigated in more de- 
tail. 

It is reassuring that our conclusion of a small spread in '^Ceff 
is reached independently by the study of isotopic s-process sig- 
natures in pre-solar silicon car bide (SiC) grains from carbon 
stars. While Eugaro et akl (l2003l) needed a spread of ~ 24 in the 
values of '■'Ceff to cover the SiC data, more recent laboratory 
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Fig. 9. Stellar population synthesis results of MS/S stars (top panels) and SC/C stars (bottom panels) compared to the observations. 
The grey scale is a logarithmic measure of the normalized number distribution of stars in the [ls/Fe]-[Fe/H] plane (left panels) and 
the [hs/Fe]-[Fe/H] plane (right panels), using '^Ceft = 4/3 (as derived in g4TT]and ^4T2l i. imin = 0.2 and /i3c,is = 0.025. The 
observational data are from the same references and with the same coding as those of Fig.|3]and Fig.|4] which have an average error 
given by the size of the cross in the upper right of each plot. 



studies (i Barzvk et ani2006l: iMarhas et alll2006h , by ehminating 
the effect of contamination of solar material, have also reduced 
the spread of the '^Ceff values needed to cover the grain data to 
exactly the same range we have found here. Although pre-solar 
grains can only give us information on the i-process taking place 
in stars with [Fe/H] > -1, such an independent check also gives 
us confidence in our results at lower metallicities. 

Now the question arises if the small spread that we find is 
compatible with the possible mechanisms currently proposed for 
the proton diffusion leading to the formation of t he cl3 pocket. 
We note that semi-conve ction dH oIIowell & lbenl ll988l) . hydro- 
dynamical overshooti ng (Herwig et all 19971) and gravity waves 
jPenissenkov & Tou3l2003) . which are three out of the four pro- 
posed mechanisms, produce a proton profile with the number of 
protons varying continuously from the envelope value to zero. 
The poss ible sh apes of the proton profile where discussed by 
iGorielv & Mowlavil (I2000i see their Fig. 10), who concluded 
that "the s-process predictions are only weakly dependent on the 
shape of the H-profile". Hence, from first principle, we would 
not expect a wide range of neutron exposures to occur in the cur- 
rent scenario. Future work should quantitatively link the shape 
of the proton profile to the range of '"^Cefr we have found in 
this study. The fourth process that has been proposed for the 



formatio n of the ^^C neu tron source is rotation ("Lang er et alj 
1999; Herwig et al.ll2003h . Also in this case the formation of 
the '^C pocket starts from a continuous proton profile, however, 
current models find that further mixing during the inter-pulse 
period coraplete ly inhibits the s-process (Herwig et al. 20031 
ISiess et aT ] l200l . Future work is need to ascertain this point. 

To explain the existence of galactic post-AGB stars that are 
.s-process enhanced at [Fe/H] ~ -1 we find that the minimum 
core mass for TDU m ust be 0.055-0 . 075 M r s maller than that 
found in the models of Karakas et alJ (|2002|) and lStancliffe et all 
(|2004|) . This is consistent with the findings of several authors 
who studied the carbon-star luminosity functions of th e SMC 
and LMC (e.g., [O roenewegen & de Jong 1993, Marig cTet al.l 
11991 llzzard et al JB004| and IStanchff^e et all 120051) and sets a 
strong constraint on future detailed evolutionary model results. 

We find that to account for the fact that all 1 1 Zr-enhanced 
post-AGB stars observed so far are also C enhanced, /Imi,, > 0.2 
in stars of mass M < 1.5 Mq. The precise value of this lower 
limit may depend somewhat on our adopted choices for mass 
loss, the AMR and the minimum core-mass for TDU. We leave 
a detailed analysis of this dependence for future study. However, 
it is clear from our results that a value of /I = cannot reproduce 
the observations, and we are therefore confident in concluding 
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Fig. 10. Population synthesis results for single-star yields compared to observations of extrinsic i-enhanced stars for different values 
of '^^Ceff. The grey scale is a normalized logarithmic distribution of [hs/ls] for a population of stars as a function of metallicity. The 
symbols indicate observed data from different stellar types and references as follows: -i- MS/S stars; x Ba II giants; O CH sub-giants; 
□ CH giants; ■ C giant; T halo CH gi ant; • halo yellow symb iotic; A halo C-ric h giant; A halo C-rich sub- giant; 7 halo N -rich dwarf, 
all previous gathered bv lBusso et alj (12001) ; O C stars from lAbia et alj (|2002|) ; ♦ lead stars gathered bv .Bisterzo et alj (2006). The 
observed data have an average error given by the size of the upper right cross in each plot. 



that stars of mass M < 1.5 M© must experience a significant 
amount of TDU. 

The amount of Zr enhancement in post-AGB stars is deter- 
mined by the amount of dredge up and the size of the '^C pocket. 
Assuming A^nin - 0.2 indicates that in our models the '■'C pocket 
must have a mass of 1/40 of that of the intershell to fit the ob- 
served data from Zr enhanced post-AGB stars. This is somewhat 
smaller than the size of the '^C pocket assumed in the models 
of G allino et al.. (1998) . which is about 1/20 of the mass of the 
intershell. 

Future detailed modelling of AGB stars will need to address 
the constraints we have set here. We are currently extending our 
analysis of extrinsic .s-enhanced stars by means of binary pop- 
ulation synthesis, i.e., by explicitly following the evolution and 
interaction processes of populations of binary stars. 

Acknowledgements. The authors are deeply indebted to Roberto GalUno for pro- 
viding the data from his detailed nucleosynthesis models and for insightful dis- 
cussions and comments on this work. ABM thanks Maarten Reyniers and Hans 
van Winckel for useful communications, discussions and for providing data on 
post-AGB stars. ABM also thanks Selma de Mink for refreshing his knowledge 
of statistical analysis. RGI is supported by NWO. ML is supported by NWO 
(VENI fellow). 



References 

Abia, C., Busso, M., Gallino, R., et al. 2001, ApJ, 559, 1117 

Abia, C., Dommguez, I., Gallino, R., et al. 2002, ApJ, 579, 817 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Aoki, W., Ryan, S. G., Norris, J. E., et al. 2002, ApJ, 580, 1149 

Barzyk, J. G., Savina, M. R., Davis, A. M., et al. 2006, submitted to Meteoritics 

& Planetary Science 
Bisterzo, S., GalHno, R., Straniero, O., et al. 2006, Memorie delta Societa 

Astronomica Italiana, 77, 985 
Boothroyd, A. I. & Sackmann, I.-J. 1988, ApJ, 328, 632 

Busso, M., Gallino, R., Lambert, D. L., Travaglio, C., & Smith, V. V. 2001, ApJ, 
557, 802 

Busso, M., Gallino, R., & Wasserburg, G. J. 1999, ARA&A, 37, 239 

CaiTaro, G., Girardi, L., Bressan, A., & Chiosi, C. 1996, A&A, 305, 849 

Cui, W. & Zhang, B. 2006, MNRAS, 368, 305 

Denissenkov, R A. & Tout, C. A. 2003, MNRAS, 340, 722 

Frost, C. A. & Lattanzio, J. C. 1996, ApJ, 473, 383 

GalHno, R., Arlandini, C., Busso, M., et al. 1998, ApJ, 497, 388 

Giridhar, S. & Arellano Ferro, A. 2005, A&A, 443, 297 

Goriely, S. & Mowlavi, N. 2000, A&A, 362, 599 

Groenewegen, M. A. T. & de Jong, T. 1993, A&A, 267, 410 

Herwig, F. 2005, ARA&A, 43, 435 

Herwig, F, Bloecker, T, Schoenberner, D., & El Fid, M. 1997, A&A, 324, L81 

Herwig, F, Langer, N., & Lugai'o, M. 2003, ApJ, 593, 1056 

HoUowell, D. & Iben, I. J. 1988, ApJ Lett., 333, L25 

Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 



Bonacic Marinovic et al.: 5 -process and population synthesis of AGB stars 



13 



Lead stars (C13gf,o=1/12) 

— I 1 r 



Two ridges of 
high probabiity 






-0.5 

-1 

-1.5 
-2 



Lead stars (C13gf,(,=1/6) 

"I 1 1 r 



-0.5 




Lead stars (C 1 3ef,= 1 /3) Lead stars (C 1 3ef,=2/3) 




-2.5 -2 -1.5 -1 -0.5 -2.5 -2 -1.5 -1 -0.5 



[Fe/H] [Fe/H] 

Fig. 12. [Pb/hs] results from our population synthesis of lead stars, shown in the same way as in fig. [TO] and compared to the data 
gathered by .Bisterzo et al., (.2006.) . The upper right cross shows the average eiTor of the data points. 



Iben, Jr., I. 1977, ApJ, 217, 788 

Iben, Jr., I. & Renzini, A. 1983, ARA&A, 21. 271 

Izzard, R. G., Dray, L. M.. Karakas, A. I., Lugaro, M., & Tout, C. A. 2006, A&A, 
460, 565 

Izzard, R. G., Tout, C. A., Karakas, A. I., & Pols, O. R. 2004, MNRAS, 350, 407 
Jorissen, A., Frayer, D. T., Johnson, H. R., Mayor, M., & Smith, V. V. 1993, 
A&A, 271,463 

Karakas, A. I., Lattanzio, J. C, & Pols. O. R. 2002, PASA, 19, 515 
Kroupa, P, Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 
Langer, N., Heger, A., Wellstein, S., & Herwig, F. 1999, A&A, 346, L37 
Lattanzio, J. C. 1989, ApJ Lett., 344, L25 

Lugaro, M., Herwig, F., Lattanzio, J. C., Gallino, R., & Straniero, O. 2003, ApJ. 
586, 1305 

Marhas, K. K., Hoppe, P., & Ott, U. 2006, Meteoritics & Planetary Science, 
Vol. 41, Supplement, Proceedings of 69th Annual Meeting of the Meteoritical 
Society, held August 6-11, 2006 in Zurich, Switzerland., p. 5248, 41, 5248 

Marigo, P, Girai'di, L., & Bressan, A. 1999, A&A, 344, 123 

Merrill, S. P W. 1952, ApJ, 116, 21 

Mowlavi, N. 1999, A&A, 344, 617 

Paczyiiski, B. 1970, AcA, 20, 47 

Pols, O. R., Schroder, K.-P, Hurley, J. R., Tout, C. A., & Eggleton, P P 1998, 

MNRAS, 298, 525 
Pont, F. & Eyer, L. 2004, MNRAS, 351, 487 

Reimers, D. 1975, Circumstellar envelopes and mass loss of red giant stars 
(Problems in stellar atmospheres and envelopes.), 229-256 

Reyniers, M., Van Winckel, H., Galhno, R., & Straniero, O. 2004, A&A, 417, 
269 

Siess, L., Goriely, S., & Langer, N. 2004, A&A, 415, 1089 
Smith, V. V. & Lambert, D. L. 1990, ApJS, 72, 387 
Stancliffe, R. J., Izzard, R. G., & Tout, C. A. 2005, MNRAS, 356, LI 
Stancliffe, R. J., Tout, C. A., & Pols, O. R. 2004, MNRAS, 352, 984 
Straniero, O., Chiefii, A., Limongi, M., et al. 1997, ApJ, 478, 332 



Straniero, O., Dommguez, I., Cristallo, R., & Gallino, R. 2003, PASA, 20, 389 

Travaglio, C., Galli, D., Gallino, R., et al. 1999, ApJ, 521, 691 

van Eck, S., Goriely, S.. Jorissen, A., & Plez, B. 2003, A&A, 404, 291 

van Winckel, H. 2003, ARA&A, 41, 391 

Van Winckel, H. & Reyniers, M. 2000, A&A, 354, 135 

Vassiliadis, E. & Wood, P R. 1993, ApJ, 413, 641 

Appendix A: Details of the synthetic TP-AGB model 

iHurley et alJ (l2000i hereafter HOO) developed a comprehensive 
synthetic evolution code for s i ngle s tellar evolution, based on 
detailed models of lPols et alJ (1 19981 hereafter P98). The P98 
models include convective core overshooting, but do not un- 
dergo thermal pulses during the AGB phase. To model TP-AGB 
stars, HOO apply the same luminosity-core mass relation as in the 
early (E-) A GB phase. The TP-AGB section of the HOO code was 
improved by llzzard et"aLl (l2004i, hereafter 104) w ith a synthetic 
code b ased on the detailed model calculations bv lKarakas et alJ 
(l2002i h ereafte r K02) which undergo thermal pulses. The latest 
version (llzzar d et al. 2006) also follows the nucle osynthesis of 
many jf-process isotopes based on the models of iGallino et alJ 
(119981) . 

However, unlike the P98 detailed models, the K02 detailed 
models do not include convective overshooting. This means that 
the transition from the E-AGB to the TP-AGB phase is not self- 
consistent in the 104 code, leading to discontinuities in the evo- 
lution of the stellar core mass, the luminosity and the radius of 
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Fig.A.l. Luminosity-core-mass relation of a 3 star with 
metallicity Z - 0.01. The arrows show the way in which the 
core evolves in time. During the E-AGB phase, the hydrogen- 
exhausted core mass remains constant at 0.716 Mq. In our model 
the core mass at the start of the TP-AGB is the same as in the 
H02 model, which is greater than in the 104 model. 
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Fig. A.2. Evolution of luminosity in time of the same star as in 
fig. lA.ll The E-AGB evolution is slightly modified in our model 
and the discontinuity in luminosity is fixed. 



Table A.l. Coeflicients aij for Eq. iAAi . 



the AGB star. To overcome these problems we have made mod- 
ifications to the 104 code, which are explained in detail below. 

A. 1 . TP-AGB initial H-depleted core mass 

For an AGB star to enter the TP-AGB phase, its helium inter- 
shell mass must be small enough to undergo a flash. This oc- 
curs at the end of the E-AGB phase when the helium burning 
shell almost reaches the hydrogen burning shell. In our code the 
H-depleted core mass at the beginning of the TP-AGB phase, 
Mc,\TP, is given by the mass of the H-depleted core at the base 
of the AGB, Mcbagb- However, stars with Mc,bagb ~ O.8M0 
and larger undergo second dredge-up by the end of the E-AGB 
phase, which makes Mc itp < Mc.bagb- 

When convective overshooting is taken into account, the 
value of McjTP for a star with a given initial mass is larger than 
that of a star with the same initial mass without overshooting. 
The effects of overshooting become important in AGB stars of 
M > 2.5Mi^ having a direct impact on the subsequent evolu- 
tion. 

We have improved the linear relation between Mc bagb and 
Mc,\TP from HOO (denoted as Mc_du in that paper) with a 
quadratic one, based on data from the P98 models (hereafter 
OV): 

Me,iTP,OV= 9.292 XlO-2Me.BAGBX(Me,BAGB + 1.983) 

(A.l) 

+ 5.865 X 10"', 

in solar units. This fit is not very good for stars with M < 1 .5Mq, 
but in these the effects of overshooting become negligible and 
the fit to the K02 models (Mcjtp,ko2) is appropriate, so we use 
instead 

Mc,np = max[ min(Mc,iTP,ov, Mcbagb), 

-. (A.2) 

A^c,lTP,K02 ■ 



' The effects of overshooting are also strong in lower-mass stars on 
the main sequence and the horizontal branc h phase, but a re largely 
wiped out during the red giant branch phase (iPols et aljl998t) . 
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Fig. lA.ll shows a luminosity-core-mass plot where the core- 
mass prescriptions of the different synthetic models can be com- 
pared. 

A.2. CO core mass and luminosity on E-AGB phase 

As in HOO, we model the luminosity during the E-AGB phase 
with a power-law relation of the form Le-agb = 'i^cc where 
Mco is the mass of the CO core during the E-AGB phase. The 
coeflicients li and I2 are found from the values of Le-agb and 
Mco at the beginning and at the end of the E-AGB phase. 

For the beginning we use Eq. (56) from HOO for the luminos- 
ity at the base of the AGB and the following relation for Mco at 
the base of the AGB 

MccBAGB = (aiM"' + aif^\ (A.3) 

which fits the data from the P98 detailed models, expressed in 
solar units. The coefficients depend on metallicity and are calcu- 
lated with the formula 

at = fl,- 1 + [Fe/H](fl,,2 + [Fe/H](fl,,3 + [Fe/H]fl,-4)), (A.4) 

where aij are the dimensionless coefficients listed in table lA.ll 
For a smooth transition to the TP-AGB phase, we require 

Mco.end = Mc,iTP and LE-AGB,end = ^ITP, (A.5) 

where Mco.end and LE-AGB,end are the CO core mass and the lu- 
minosity at the end of the E-AGB phase, respectively. Lup is 
the luminosity at the beginning of the TP-AGB phase and is dis- 
cussed in ^A.3I A comparison between the time evolution of the 
luminosity described by the different synthetic models is shown 
in fig.O 



Bonacic Marinovic et al.: 5 -process and population synthesis of AGB stars 



15 



4.1 
4 
3.9 
3.8 
3.7 
3.6 
3.5 
3.4 
3.3 
3.2 



1 

- 


1 


1 1 


1 


1 1 












I 




t. 


"7 
1 


/ 
/ 
f 


i 
I 
I 


1 

1 




1 


f 
1 










1 

— \ 

1 


1 


1 1 


M=1.0 Mq 
Z=0.004 

1 1 1 





0.56 0.565 0.57 0.575 0.53 0.585 0.59 0.595 0.6 
"cNODUP [l^oi 

4.2 



4 
3.8 
3.6 
3.4 
3.2 



T 1 r 




J L 



0.54 0.55 0.56 0.57 0.58 0.59 0.6 0.61 




2.8 

0.54 0.545 0.55 0.555 0.56 0.565 



4.25 
4.2 

4.15 
4.1 

4.05 
4 

3.95 
3.9 

3.85 
3.8 




M=3.0 Mp 
Z=0.004 



0.7 0.71 0.72 0.73 0.74 0.75 0.76 0.77 0.78 0.79 
Mc.NODUP [Mq] 



1 — I — I — I — I — I — I — r 




0.69 0.7 0.71 0.72 



M=3.0 M„ 
Z=0.02 



J I I I I I I I L 



3.2 

0.56 0.58 0.6 0.62 0.64 0.66 0.68 0.7 0.72 0.74 0.76 
Mn.NOniJP [MJ 




M=5.a Mp 

z=a.oo4 



0.9 0.905 0.91 0.915 0.92 0.925 0.93 0.935 0.94 



™c,NODUP l">el 




0.87 0.88 0.89 0.9 0.91 0.92 0.93 
l^cNODUP [MqI 

— r 



- =i 4.4 - 




0.85 0.9 0.95 1 1.05 



Fig. A .3. Examples of our analytic fit (solid lines) compared to the data from the detailed models (dashed Unes) from lStancliffe et alJ 
(12004 for different stellar masses and metallicities, as indicated in the panels. 



A.3. Luminosity during TP-AGB phase 



up: 



Following 'P aczynskil (Il970h . it has become usual practice to 
apply a linear luminosity-core-mass relation in synthetic AGB 
models. However, this relation fails in the presence of dredge-up 
and hot bottom burning. This becomes a problem when com- 
paring models which include overshooting to models that do 
not: two TP-AGB stars with the same core mass, but different 
envelope mass, do not behave in the same way. To overcome 
this we design a luminosity relation which depends on both the 
core mass and the envelope mass of the star It successfully re- 
produces both the P98 models which include overshooting and 
models calculated with the same code used by P98 which do not 
include overshooting. Stancliffe et al.l (1200 4, hereafter S04) em- 
ployed a modified version of the code used by P98 and computed 
detailed models which do not include overshooting, but do un- 
dergo thermal pulses. In an attempt to make our rapid TP-AGB 
models more self-consistent we have modified our luminosity 
relation to fit the S04 models, assuming that with the inclusion 
of overshooting the S04 models are modified in the same way 
as those of P98. In addition to this, we have found that in the 
presence of dredge-up the results of the detailed stellar evolu- 
tion models are better represented by relating the luminosity to 
^c.NODUP, the core mass as it would be in the absence of dredge- 



logio(ic,NODUp) = ( logio(15300 + 9440Z) 

3.45 logio(Mc,NODUP + 0.1) 



3Mc,N0DUP - A^clTP 
0.045 

+ 0.08, 



(A.6) 



A^cNODUP / 



in solar units and where Z is the metallicity. In stars with a mas- 
sive convective envelope (Mgnv S 2.5 M©), the temperature at 
the bottom of the envelope is high enough for hot bottom burn- 
ing (HBB), which enhances the luminosity of the star by a factor 
/hbb due to the enhanced hydrogen abundance in the H-burning 
shell. The expression 

Me„v,iTp(l - 0.53 [Fe/H]) 



logio(/HBB) = max < 1.3 log 



- 0.08 , 



env.HBB 



(A.7) 

fits /hbb X Lc.NODUP to the luminosity of the S04 models, where 
Menv.iTP is the mass of the envelope at the beginning of the TP- 
AGB and we have set the threshold mass for the onset of HBB 

Menv.HBB = 2.4M0. 
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The S04 models have no mass loss, so we have to use a differ- 
ent set of models to fit the behaviour of the luminosity in terms 
of the loss of the envelope mass. Based on the K02 models the 
contribution of the envelope to the luminosity can be modelled 
as 



logio(/env) = max { logio(/HBB) + 



0.46 min 



logio 



Me, 



Me, 



,0 



, 



(A.8) 



The expression /env x Lc.nodup matches the asymptotic be- 
haviour of the luminosity observed in the S04 detailed models, 
which is reached after several pulses. To model the first few 
pulses we write the luminosity as 



L = 



-^c.NODUP 



X /ei 



-^TPl 



/ 



(A.9) 



where Ltpi is fitted to the S04 models with the expression 



logio(LTPl) 

and 



-1 .8427m2.^pj -h5.3861Mc,tpi -hO.97823 (A. 10) 



/= 1- 



1 - 



M, 



c.NODUP 



M< 



c.lTP 



0.075 



25/3 



(A.ll) 



accounts for the luminosity growth during the early pulses to its 
asymptotic value Lc.nodup x /env In Eqs. ( lA.lOl l and (lA.lll i all 
quantities are expressed in solar units. Fig. lA.3l shows a compar- 
ison between our fit and the S04 detailed models. In the worst 
cases, our fits have an error smaller than 10% in luminosity, 
which is less than the difference between the S04 and K02 de- 
tailed models. 



A.4. Radius during TP -AGB pinase 

We use the HOO radius-luminosity relation, given that it accu- 
rately fits the S04 models, but we have slightly modified the in- 
terpolation procedure. The HOO radius-luminosity relation (Eq. 
(74) in that paper) was designed to fit the P98 models at seven 
different metallicities in the following way: 

/?AGB = ACL*' + Z72i'0, (A.12) 

where L is the luminosity, A depends on the mass, and bi, on 
the metallicity. For intermediate metallicities the radius was cal- 
culated after interpolating linearly in Z between the coefficients 
that fit these seven metallicities. Some of these coefficients are 
exponents and we found the results for interpolated metallici- 
ties to behave erratically. We therefore interpolate linearly in Z 
between the radius values calculated with Eq. lA.12l for the adja- 
cent known metallicities, which is in much better agreement with 
detailed models at intermediate metallicities calculated with the 
same code as used by P98. 



